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Abstract 

In this paper, we study the dynamics of a hohow spherical matter cohapsing with very large 
initial velocity. The spacetime is initially very similar to the Vaidya solution, and the deviations 
from this background are treated perturbatively. The equations of state for radial pressure = 
kp and tangential one = wp with constant k and w are assumed. We find for the case of 
equations of state k < \ and < u; < 1 that the initial velocity, which is nearly the speed of 
light, is strongly decelerated. This result implies that the pressure is essential to the property 
of singularity formation in gravitational collapse even for initially nearly light-speed collapse. By 
contrast, in cases with the negative tangential pressure, the present result implies that the central 
naked singularity similar to that of the Vaidya spacetime can be formed, even though the radial 
pressure is positive, and the weak, strong and dominant energy conditions hold. Especially, in the 
case of u; < —(1 — A;)/4, the high-speed collapse will produce the spacetime structure very similar 
to that of the Vaidya spacetime. 
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I. INTRODUCTION 



Einstein presented the theory of general relativity which describes the gravitational force 
in terms of the spacetime curvature. He formulated the field equations which relate the 
geometry of the spacetime to matter fields. The earliest exact solutions of these equations 
were the Schwarzschild metric representing the exterior of a spherically symmetric star and 
the Friedmann cosmological models. These solutions have spacetime singularities where the 
energy density or spacetime curvature diverges and the usual description of the spacetime 
is impossible there [1]. 

It is a well known phenomenon that gravitational collapse of massive objects results in 
the formation of spacetime singularities in our universe. There are two kinds of spacetime 
singularities from a point of view of their visibilities. A spacetime singularity is said to be 
naked when it is observable to local or distant observers, and the remaining one is said to 
be covered. The singularity theorems of Hawking and Penrose show that the formation of 
spacetime singularities is not rare in our universe if the general relativity is correct and the 
matter or radiation fields satisfy physically reasonable energy conditions j^. However, these 
theorems do not provide information about the visibility of the spacetime singularity. 

About the visibility of the singularity, Penrose proposed a conjecture called the cosmic 
censorship conjecture which has two versionsjsl. The weak version states that the spacetime 
singularities produced by gravitational collapse of physically reasonable matter fields, which 
develops from generic non-singular initial data, are always covered by horizon, whereas the 
strong version claims that there is no singularity visible to any observers. Some rather 
serious counterexamples have been found for the strong version 0, 0]. As for the weak 
version, any precise theoretical or mathematical proof has not yet been given, although it 
has many physical applications in black hole and other areas in astrophysics. This motivates 
that a detailed study of dynamically developing gravitational collapse models is necessary 
to obtain a correct form of the cosmic censorship. Several examples (0]-[l4| and references 
therein) have been studied so far which admit both black hole and naked singularity solutions 
depending on the choice of the initial data. Most of the work on gravitational collapse has 
been done by considering dust fluid due to the existence of an exact solution. However, the 
assumption of dust fluid might be too restricted as the effects of pressure can not always be 
neglected in the formation processes of the spacetime singularities. Thus it is important to 
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discuss this issue by including pressure. 

Ori and Piran l^-[l7| investigated self-similar spherically symmetric perfect fluid collapse 
by assuming the equation of state p = kp. They found that a naked singularity is formed 
for < A; < 0.0105. They have also shown that there exist naked- singular solutions with 
oscillations in the velocity field for < k < 0.4. Later, these results were extended for 
< k < 0.5625 by Foglizzo and Henriksen ilSjl. The same results were also provided without 

□ n 

self-similarity assumption by one of the present author TH and Maeda[19|-[21|]. Giambo et al. 
have investigated naked singularity formation in perfect fluid collapse without self-similarity 
assumption analytically 221]. Goswami and Joshi 23|] have also investigated analytically the 
local geometry near the central shell focusing singularity formed by a spherical collapse of 
a perfect fluid with the equation of state in the form p = kp and have shown that the 
initial condition is crucial for whether this singularity is naked. However, it is difficult to 
construct a general global solution analytically. Towards the progress of analytical studies, 
it is important to develop a new approximation scheme (analytical procedure) to discuss 
gravitational collapse which leads more definite results. 

One of the present authors KN and Morisawa studied the cylindrically symmetric gravi- 
tational collapse of a thick shell composed of dust by introducing high-speed approximation 
scheme 24] . The same authors generalized this work for the perfect fiuid case 25|] . In a recent 
paper 26| . two of the present authors MS and ZA have extended this work by considering two 
perfect fluids. These investigations have provided interesting results about the gravitational 
collapse. It would be worthwhile to explore whether these results hold for spherical collapse 
or not. This motivated us to develop a high-speed approximation scheme for spherically 
symmetric system with the type I matter 27|. There are many studies about the naked sin- 
gularity formation by the same type of matter in the spherically symmetric system. Dwivedi 
and Joshi showed by the local analysis that the initial data is crucial for whether the naked 
singularity forms at the symmetric center by the gravitational collapse of general type I 



matter 
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and Joshi 



_ lecently, this work extended to the higher dimensional spacetime by Goswami 



291 ]. The spherically symmetric matter with vanishing radia! 
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pressure and non 
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33, 



35]. 



vanishing tangential pressure has been studied by various authors 
In this paper, we consider the case treated by Ref.[28i] but we construct global analytic 
solutions by using high-speed approximation. 

The paper is organized as follows. In section II, we write down the Einstein equations 
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for the spherically symmetric spacetime with a type I matter in the single null coordinate 
system. The null dust solution is investigated in section III. Section IV is devoted to discuss 
the high-speed approximation scheme for the general type I mater. The effects of pressure 
on the high-speed gravitational collapse are discussed in section V. Finally, the summary of 
the results is given in section VI. 

In this paper, we adopt the geometrized unit, i.e., c = 1 = G and follow the convention 
of the Riemann and metric tensors and the abstract index notation adopted in the textbook 



by Wald [36|; the latin indices denote the type of a tensor, whereas the Greek indices denote 



the components of a tensor. 

II. SPHERICALLY SYMMETRIC MATTER WITH ANISOTROPIC PRESSURE 

We focus on the spacetime with spherical symmetry. For later convenience, we adopt the 
single null coordinate system in which the line element is given by 

ds^ = -A{v, r)dv^ + 2B{v, r)dvdr + {d9^ + sin^ Odcp^) , (1) 

where we assume that B is positive, and this assumption implies that the coordinate v is 
the advanced time, i.e., a radial curve of constant v is future directed ingoing null. Then 
the non-trivial components of the Einstein equations are given by 

= S7iT\, (2) 
= SttT^, (3) 
= SnT\, (4) 

= SvrT^, (5) 
1 / A"\ B' f ■ A'\ A /A' B'\ , 




where the dot denotes the derivative with respect to the advanced time v, while the prime 
denotes the derivative with respect to the radial coordinate r. 

We study the dynamics of the spherically symmetric type I matter whose stress-energy 
tensor is 2?! 

T\ = pu^Ub + PRS^Sb + PTn\, (7) 



where u°' is 4-velocity of a constituent particle, s° is the unit radial vector normal to and 
Qf^h is defined by 

n\ = 6\ + u'^Ub - s'^Sb, (8) 

and thus p, pn and pt are the energy density, radial pressure and the tangential pressure, 
respectively. 

We write the components of and Sa in the forms 

u'' = N {V, -1 + V, 0, 0) , (9) 
= NB {1-V,V,0,0), (10) 



where 



N= ^ (11) 
y^V {2B + V{A-2B)} 



We define new variables D, Pr and Pt as 



^ N^V^{p + Pr) ^ r^B{p + pn) 

sine V {2B + V{A-2B)}' ^ ' 



^' sin^ V{2B + V{A-2B)}' ^ ^ 

.= N'V^P^ = !^ (14) 

^ ■ sin^ V{2B + V{A-2B)}' ^ ' 
where g is the determinant of the metric tensor. Then the stress-energy tensor is written in 
the form 

T\ = [Dk^'h + V{2B + V{A - 2B)} {Pr5^ + (Pt - Pr)^^}] , (15) 
where the components of the vector field k"- are 

k' = j^ = iy^-^ + V,0,0). (16) 



III. NULL DUST LIMIT 
A. Metric 

It is easy to see that in the limit of — > with D, Pr and Pr fixed, the stress-energy 
tensor (|T5l) becomes that of the null dust which belongs to the type II matter j27||. 
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where in this hmit we have 

P = (0,-1,0,0). (18) 
We can easily check that B~^k^ is the tangent of the ingoing null geodesic, i.e., 

k'^Va [B~^h) = 0. (19) 

Thus the equation of motion V aT°'h = leads 

Va {r-'^Dk") = 0. (20) 

The above equation reduces to 

(BD)' = 0, (21) 

and thus we find that BD is a function of the advanced time v only. Further since Eq.([3]) 
reduces to B' = 0, we have B = B{v). Therefore, we have D = D{v). If we introduce a new 
advanced time v defined by 

v = j B{v)dv, (22) 

we have a new line element 

ds^ = -A{v, r)dv^ + 2dvdr + {d9^ + sin^ Odcp^) , (23) 

where A = A/B"^. Thus without loss of generality, we can assume B = 1 and will do so 
below in this section. 

The remaining non-trivial components of the Einstein equations become 

^-i(l-A) = 0, (24) 

A ^ (25) 

r 

A" A' 

— + - = 0. (26) 

From Eq.( l2^ . we obtain 

^ = l-?M!i. (27) 

r 



From Eq. (!25!) and the above equation, we have 

dM 



47rD{v). (28) 



dv 

If we fix the dependence of D on v, the solution is completely determined. This solution is 
known as the Vaidya solution. 



B. Central singularity 



The precursory singularity in spherical gravitational collapse formed at the symmetry 



center can be naked singularity 



37( 1 ■ The central singularity in Vaidya spacetime can be 



naked if the mass function M{v) satisfies some condition 38|, |39|, |4l|. Here we give a brief 



review about this issue by following the analysis in Ref. [39|. 

The future directed ingoing null condition is v =const, whereas the future directed out- 
going null condition is given by 

dr 1 ^ 2M\ , , 

1 - — • (29) 



dv 2 \ r ^ 

We assume that -D(f) has a compact support so that we have M[v) = for f < and 
M{v) 7^ for f > 0. This assumption guarantees that the symmetry center r = is regular 
initially, i.e., for v < 0. The central singularity forms at f = = r. To know whether the 
central singularity is naked, we investigate the existence of the outgoing null geodesies from 
V = = r. 

Suppose that the r = ro{v) is a solution for Eq.( l29i) which emanates from the central 
singularity v = = r. The mass function is written by using this solution as 

Miv) = ^roiv) [1 - 2roiv)] . (30) 

The above equation is just the condition on the mass function for which the central singu- 
larity becomes naked. It should be noted that there is only one ingoing null geodesic v = 
that hits the central singularity v = = r. Hence, if there is a one-parameter family of null 
geodesies emanating from the central singularity, the central naked singularity is null, while, 
if there is only one, the central singularity is instantaneous, i.e., an event in conformally 
extended spacetime manifold. Solutions in such a one-parameter family are written in the 
form 

r{v) = ro{v) + z{v;a), (31) 

where a parametrizes the solutions. Substituting the above form into Eq.( !29l) . we have 

dz ^ z{2ro - 1) 

dv 2{z + ro) ■ ^ ^ 

We search for the solutions which behave as 2; ^ for v 0. 

Suppose ro ~ Pv" near the central singularity, where both of a and /3 are positive. Then 
we find following three cases: 



FIG. 1: The conformal diagram of the Vaidya spacetime in the cases (i) and (ii). The central naked 
singularity is null and thus there is a family of the future directed outgoing radial null geodesies. 



(i) M(f) ~ /3t;"/2 for a > 1; 

(ii) M{v) ~ 13(1/2 - (3)v for a = 1; 
(ill) M{v) for a < 1. 

For case (i), the solution for Eq. (!32l) is given by 



z av "exp ( — r I and 2; ~ -. (33) 

^ ' 2(3{a - l)i;"-i / 2 ^ ^ 



For case (ii), we have, for /3 7^ 1/4, 

J ir\7.\ 1-2/3 -I- 

1-4/3 



, I 2/3 2z , 
V ~ a\z\ 1-2/3 + — , (34) 



while we have, for f3 = 1/4, 

V ^ az + 4zln\z\. (35) 

We see from Eq.( !34l) that (3 should be less than 1/2 so that 2 for f 0. This condition 
guarantees that the mass function M{v) is non- negative near v = 0. For case (iii), the 
mass function M{v) is negative near f = 0. Thus in this case, the spacetime singularity of 
f > and r = is necessarily naked but it is non-trivial whether the central singularity at 
f = = r is null. We can see that the solution for Eg. (1321) is only z ~ —2(3v°', and from 
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FIG. 2: The conformal diagram of the Vaidya spacetime in the case (iii). The central naked 
singularity is instantaneous. The final product of the gravitational collapse is the Schwarzschild 
spacetime with negative mass. 




FIG. 3: The conformal diagram of the Vaidya spacetime in the case that the central singularity is 
not naked. The mass function M behaves as M ~ fiv'^ with /i > 1/16 and 7=1 near v = 0. 
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Eq. (l3T!) . we have 

r ~ < 0. (36) 

Since r should be non-negative, this is an unacceptable solution. Thus for case (iii), the 
solution r = ro(f) is the unique outgoing radial null geodesic which emanates from the 
central singularity. Further the radial curve f = is the unique future directed ingoing radial 
null which hits the central singularity. Thus we may conclude that the central singularity 
of negative mass function is instantaneous. 

We see from the above results that if the mass function M behaves as M ~ iiv'^ near 
f = with /i > 1/16 and 7=1, then the central singularity is not naked. 

We depict the conformal diagrams for the cases (i) and (ii) in Fig.l, while for the case 
(iii) in Fig.2. The covered case is depicted in Fig. 3. 



IV. HIGH-SPEED APPROXIMATION SCHEME 

The high-speed approximation is a kind of linear perturbation scheme in which V is a. 
small variable of order e. Then we write 

A=l-^ + M^,r), (37) 
B = l + 5B{v,r), (38) 
D = DB{v)[l + 5D{v,r)l (39) 

where 5a, Sb and 6d are also assumed to be 0(e), and 

Mb = ^nDB. (40) 

We take into account the terms up to the first order of e. Hereafter, we assume that V is 
non-negative so that is a causal vector. As in the previous section, we also assume here 
that Db{v) has a compact support so that we have Mb{v) = for f < and Mb{v) ^ for 
v>0. 

From Eq. (fT5|) . we have 

T\ = (41) 
This equation and Eq.Q lead B' = 0{e^), and thus we have 

B = l + 5B{v) + 0{e^). (42) 
10 



The above result means that B is the function of only v up to the first order of e, and hence, 
without loss of generality, we can assume i? = 1 up to this order by the same reason as in 
the case of the Vaidya solution. Then Eqs.([2]) and agree with each other up to the first 
order of e. 

The remaining apparently independent components of the Einstein equations of 0(e) are 



Ha)" 



-8n{DB - 2Pn)V, 



327T 



-PtV. 



The equations of motion for the matter VaT^^ = of order 0(e) are given by 



[{DB-2PK)Vy+DB 



0, 



(43) 
(44) 

(45) 



(46) 
(47) 



[{Db-2Pk)V]' + -PtV = 0. 

Eq. (l45!) is not necessary, since this equation is derived from Eqs.f l43p and (147!) . Differentiating 
Eq. fHSl) with respect to r, we have 

{r6A)" = -87r[{DB-2Pn)V]'. (48) 

Substituting Eq. (l471) into the right hand side of the above equation, we have Eq. (l45l) . Eq. (H6|) 
is also not necessary, since this equation is derived from Eqs.( l43|) and f j44l) . Differentiating 
Eq. (H3l) with respect to v, we have 

(rSA)' = -8n[{DB-2P^)Vy . (49) 

Differentiating Eq.( l44l) with respect to r, we have 

{r5Ay = SttDb 



(50) 



1 + 



From Eqs. fHOj) and fl50l) . we have Eq. fj46|) . Thus independent equations are only Eqs.fHS 
(HID and ^ 



V. THE ANALYSIS OF THE HIGH-SPEED COLLAPSE 

In order to write down formal solutions for the perturbation variables, we introduce new 
variables k{v,r) and w{v,r), defined by 

p^_ = k{v,r)p and = w{v,r)p. (51) 
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Here, we briefly review the so-called weak, strong and dominant energy conditions, which 
might be satisfled by physically reasonable matter. These conditions are rewritten in the 
form of the conditions on the values of k and w. The weak energy condition leads the 
following conditions on p, and p^, 

p > 0, p + Pr>0 and p + PT>0; (52) 

the strong energy condition adds one more condition 

p + PR + 2pT>0; (53) 

the dominant energy condition leads further two conditions, 

P > \pr\ and p > \pt\- (54) 

Assuming that p is non-negative, all the three energy conditions are guaranteed, if and only 
if following conditions for k and w are satisfled: 

-1<A;<1, -l<w <1 and l + k + 2w>0. (55) 

Hereafter we assume the non-negativity of p. It is worthy to notice that both of k and w 
are bounded above and below by virtue of the energy conditions. 
From Eq. llST]) . we have 

Pr PR k Pt Pt w 

-TT = = 1 7 and — = = -. (56) 

D p + pr k+1 D p + p^ k + 1 ^ ^ 

Using the above equation, Eq. (l47I) becomes 

'~'v\ + ^V = 0, (57) 



(58) 



1 + k J r{l + k) 
where we have used the fact Db = Dsiy). From the above equation, we have 

4w 



In 



l + k 



r(l - k)' 



Clearly, the high-speed approximation scheme is not applicable to the case oi k = ±1: the 
upper sign corresponds to the stiff matter and the lower sign corresponds to the cosmological- 
constant-like matter. Hereafter we assume k ^ ±1. 
The formal solution of Eq. (1581) is given by 

F = C(i;)exp -/ TT^— ' (59) 



1 + /c \ J \ — k{v, x) X 
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where C{v) is an arbitrary function of v but it is set so that V is positive at least initially. 
Substituting the above equation into Eq. fH3l) and integrating it, we have 



exp 



1 — k{v, x) X 

where we set the integration constant so that 5a is finite in the limit of r 
Then substituting Eqs. (j59l) and (!60l) into Eq. (l44l) . we have 



(60) 
if possible. 




+ 



X exp 



Aw{y,x) dx 
1 — k{v^ x) X 

4w(t>, z)k{v, z 
\-k{v,z) ' [1 - k{y,z)f 




4ti;(t>,x) dx 



(61) 



1 — k{v^ x) X 

Once k and w are determined, we can know the behavior of the first order perturbations by 
performing the integrations in the formal solutions (l59i) . ( l60l) and ( 16T]) . 

The behavior of the solutions fl59|) - fl6T|) near the spacetime singularity will be determined 
by the asymptotic values of k{v,r) and w{v,r) in the limit that the hollow sphere shrinks 
to its symmetry center. Thus the solutions obtained by assuming the constancy of k and w 
will give us sufficient information about what we would like to know. Hereafter, we assume 
that k and w are constant. Integration in Eq.( l59l) is then easily performed and we obtain 

1 + k 



V 



1 - k 



C{v)r 



4:W 

' 1-k 



(62) 



Using the above result, integrations in Eqs. (l5Ul) and (1^ are also easily performed, and we 
obtain 



-8nCDBF{r;k), 



DbSd = CDb ( + ^ ) r-i^ + {DbC + DBC)F{r- k), 



where 



F{r;k) 



(1 - k - 4:w)'^{l - k)r i-fe for 1-k-Aw^O, 
r~^lnr for l — k — Aw = 0. 



(63) 
(64) 

(65) 



From the causality requirement 401], we assume k < 1 and w < 1. Then we study the 
following three cases, < w < 1, — (1 — A;)/4 < w < and w < — (1 — fc)/4, separately, 
below. 
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A. <w <1 

It can be seen from Eq. fl62|) that V diverges in the hmit of r 0. This imphes that, in this 
case, the high-speed collapse is necessarily decelerated by the pressure effect so significantly 
that the high-speed approximation breaks down before the singularity formation. This is a 
somewhat unexpected result, since, at first glance, one might infer that the pressure effect 
would be negligible in the situation with the strong gravity near the spacetime singularity. 



B. -(1 - A:)/4 < < 

It is easily seen from Eqs. (1621) and (l63l) that both V and Sa are finite in the limit of 
r ^ 0. Eq.§^ leads 

2Mfi „ iw 

DbSd -CDbt"^ near r = 0. (66 

r 

Here note that by assumption, we have 

0,-J^<l. (67, 

Mb/t diverges in the limit of r — with v fixed in the domain with non- vanishing Mb- 
Thus, due to Eg. (1671). the density perturbation DbSd also diverges in the limit of r — 
within this domain, if it initially does not vanish there. The high-speed approximation in 
this domain breaks down when the matter particles approach the symmetry center r = 0. 

Since Mb vanishes at f = 0, the behavior of the density perturbation DbSd in the 
neighborhood of the central singularity might be different from the above case. The behavior 
of Mb/t near the central singularity has already been shown in Sec. III. Our special interests 
are in the cases (i) and (ii), since the the background central singularity has the extent in 
the future directed ingoing null direction, while in the case (iii) and in the covered case, the 
central singularity is instantaneous (see Figs. 1-3). Thus we focus on the cases (i) and (ii) 
in which the mass function Mb behaves as Mb ~ /if" with a > 1. In these cases, Mb/t 
is finite at the central singularity in the background Vaidya spacetime. Thus the central 
singularity will be very similar to the central naked singularity of the Vaidya solution. 

Here it should be noted that this case includes a case of dust k = w = The dust is 



described by the Lemaitre-To^ 
can also be naked 42, 43, 44 



man-Bondi solution and it is well known the central singularity 



45| . The present result implies that the central singularity 
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formed by the gravitational collapse of a hollow dust sphere is very similar to that formed 
by the null dust, whereas the non-central singularity is not so. The latter is in contrast to 
the cylindrically symmetric case in which the cylindrical hollow dust collapsing with very 



high speed is well described by the null dust solution even at the spacetime singularity 



24]. 



C. w<-{l-k)/A 

In contrast to the case B, the perturbation variables V and 5a vanish in the limit of r ^ 
with V fixed. In the case oi w < —(1 — fc)/4, DbSd also vanishes in the same limit, whereas 
it is finite in the case of w = —(1 — A;)/4. Thus the spacetime structure in the neighborhood 
of the singularity at r = is very similar to the Vaidya solution. Even if the radial pressure 
is positive, the high-speed approximation is consistent until the spacetime singularity forms. 
The consistency of the high-speed approximation depends not on the radial pressure but on 
the tangential one. Here note that all of the energy conditions Eg. (1551) hold only if 



holds. 



k>-l (68) 



Recent observations imply the acceleration of cosmic volume expansion |46l. 1471 . |48|, and 
this means the existence of unknown matter components with violation of the strong energy 
condition. Thus it might be important to consider cases with the violation of the energy 
conditions as special examples of case C. A case of k smaller than —1/3 corresponds to the 
so-called dark energy, and thus the result obtained here implies that the spacetime singularity 
formed by spherically symmetric high-speed collapse of the dark energy is similar to that of 
the Vaidya solution. Further, it is worthy to note that in the phantom energy case k < —1 



which is the special case of the dark energy 49|], -Ds is negative and thus the mass function 



Mb is also negative by Eq.( l40l) . Hence the first order solution of < — 1 implies that the 
high-speed collapse of the phantom energy forms the timelike singularity similar to that in 
the Schwarzschild spacetime with a negative mass (see Fig. 2). 

VI. SUMMARY AND DISCUSSION 

Gravitational collapse is one of the most important topics in gravitational physics. The 
cosmic censorship conjecture provides major motivation to study this issue. Since there is 
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no theorem proving or disproving this conjecture or no theorem stating the generic feature 
of physical spacetime singularities, it is interesting to investigate this issue in the situation 
different from previously studied ones. 

This paper continues to study this issue and provides an extension of the previous work on 
high-speed cylindrical collapse of perfect fluid [25] to spherically symmetric spacetime with 
type I matter. To see the pressure effects on high-speed approximation scheme, assuming 
that the energy density p is non-negative, we have studied a linear equation of state for the 
radial and tangential pressures, i.e., = kp and = wp with constant k and w, in detail. 
By the causality requirement, we have restricted our attention to the case of A; < 1 and 
w <1. (The causality requirement implies k <1. The reason why k = 1 is excluded is that 
the high-speed approximation is not applicable to the cases of k = ±1.) 

One might think that the collapsing speed becomes very large due to the strong gravity 
just before the formation of the spacetime singularity. However, we have found that, in 
the case of positive tangential pressure < w < 1, the large initial imploding velocity is 
necessarily decelerated by the pressure effect and thus the high-speed approximation scheme 
becomes invalid before the singularity formation. 

In the case of the negative or vanishing tangential pressure w < 0, the behaviors of the 
perturbation variables are different from the case of < w < 1. In the case of —(1 — 
k)/A < w < 0, all of the perturbation are finite at the central singularity, if the background 
central singularity is null and naked, although the density perturbation blows up at the 
non-central singularity of the background Vaidya spacetime. This result implies that the 
central singularity will be null and naked like as that of the Vaidya spacetime, for the case 
of —(1 — k)/A < w < 0. In the cases of w < —(1 — k)/A, all of the perturbation variables are 
finite everywhere. This result implies that, in this case, the geometrical structure near not 
only the central singularity but also the non-central singularity is very similar to that of the 
Vaidya solution. Here it is worthy to note that as long as k > —1/3, this case satisfies all 
of the physically reasonable energy conditions, i.e., the weak, strong and dominant energy 
conditions. The result obtained here strongly suggests that the spacetime structure reahzed 
by the spherical matter with large enough tangential tension is well described by the Vaidya 
solution even if the physically reasonable energy conditions are satisfied. 

It is a remarkable result that the consistency of the high-speed approximation depends 
not on the radial pressure pr but on the tangential one pt- In the case of < 0, the 
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gradient of radial pressure will not stop the high-speed collapse. Thus it is reasonable that 
the consistency of the high-speed approximation depends on the only tangential pressure. 
It is non-trivial that even if A; > 0, the consistency of the high-speed approximation also 
depends on the only tangential pressure. At first glance, the gradient of the radial pressure 
with A; > seems to affect the gravitational collapse, but this is not true. The reason is 
that if the sound speed in the radial direction is less than the speed of light {0 < k < 1), 
the spacetime singularity formation by the nearly light-speed collapse can be completed 
before the effect of radial pressure gradient spreads out, since the sound cone is significantly 
narrowed down in the frame in which the matter moves with nearly the speed of light. This 
is also the reason why the high-speed approximation is not applicable to the case of k = 1. 
In this case, the sound speed in the radial direction is equal to the speed of light and thus the 
sound cone in the radial direction is equivalent to the light cone which is Lorentz invariant. 
As a result, the effect of the radial pressure gradient can spread out before the singularity 
formation by the nearly light-speed collapse is completed. The perturbative construction 
of approximate solutions on the background null dust solution is impossible in the case of 
k = l. 

The dark energy case, k < —1/3 and w < —(1 — k) /4, might be important in connection 



to the issue of the accelerated cosmic expansion |46l. l47l . |48|, although the strong energy is 
not satisfied in this case. These are included in the case of w < —(1 — k)/A. The repulsive 
gravity due to the dark energy cannot decelerate the high-speed collapse, and the formed 
central singularity can be null, and the non-central part of the spacetime singularity will be 
spacelike. As mentioned in the above, the pressure gradient force of the dark energy will 
also not stop the collapse but rather accelerate it. This property of the dark energy will be 
the reason why the high-speed collapse becomes a good approximation at the singularity 
formation in the present case and also why black holes may grow self-similarly due to the 
accretion of dark energy in the accelerated universe 50|, l51 |. 

Finally, it should be noted that the present results are valid up to the only first order 
and could be modified by the higher order effects. Thus we need to investigate the higher 
order, but this is a future work. 
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